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Abstract. Some analyzes of the apparent magnitude-redshift data of type la supernovas indicate that the suspected dark 
energy in the universe cannot be regarded as a cosmological constant of general relativistic origin or as the vacuum energy 
encountered in quantum field theories. If this is the case, our knowledge of the physical world remains deficient since no tested 
theory involves such a dark energy. Under this circumstance, an equation of state of the form p = wp is not well-motivated and 
one is unable to use the Einstein equation in this case as well. I argue that the very method of analysing the data by assuming 
exotic energy densities with strange equations of state itself is misleading and the reasonable remaining option is to make 
a model-independent analysis of SNe data, without reference to the energy densities. In this basically kinematic approach, 
we limit ourselves to the observationally justifiable assumptions of homogeneity and isotropy, i.e., to the assumption that 
the universe has a RW metric. This cosmographic approach is historically the original one in cosmology. The analysis was 
performed by expanding the scale factor into a fifth-order polynomial, an assumption that can be further generalized to any 
order. The values obtained for the present expansion rates h, q$, ro etc. are relevant, since any cosmological solution would 
ultimately need to explain them. 

Using this method, we address an important question relevant to cosmology: Was there a decelerating past for the universe? 
To answer this, the Bayes's probability theory is employed, which is the most appropriate tool for quantifying our knowledge 
when it changes through the acquisition of new data. The cosmographic approach helps to sort out models which were 
always accelerating from those which decelerated for at least some time in the period of interest. Bayesian model comparison 
technique is used to discriminate these rival hypotheses with the aid of recent releases of supernova data. It is argued that the 
lessons learned using Bayesian theory are extremely valuable to avoid frequent U-turns in cosmology. 



Modern cosmology began with the Hubble's discovery that there is a systematic relationship between the apparent 
magnitudes and redshifts of galaxies. He observed that as galaxies become fainter and fainter, their redshifts increase 
systematically. When these redshifts are interpreted as originating from a Doppler-like effect, one concludes that 
farther galaxies recede from us with greater speeds. Assuming that we do not enjoy any preferred position in the 
universe, this discovery naturally leads to the concept of an expanding universe. Ever since 1929, the expansion is the 
primary piece of information in cosmology. In this talk, my starting premises are that of an expanding, homogeneous 
and isotropic universe. 

The Friedman solution of Einstein equation in cosmology, which assumes that at large scales the universe has a 
Robertson- Walker metric and is filled with ordinary relativistic or dust-like matter, points to an instant of creation 
or initial singularity. This is usually referred to as the big bang origin of the universe. The Hubble's discovery was 
in fact an observational verification of the expanding universe model. Later observations, such as the abundance of 
light elements in the universe, the cosmic microwave background radiation (CMBR), etc. were considered to endorse 
this picture. However, this idea has found many ups and downs and many U-turns in the past eight decades and is an 
example of the curious way scientific theories evolve. 

A lacuna in modern cosmology is that the errors in the measurement of the present values of the expansion rates, 
related to the Hubble parameter Hq, deceleration parameter q$, etc., of the universe are notoriously high. For example, 
Hubble's own measurement of Hq is higher than the presently favored values by an order of magnitude. Even the 
most optimistic measurements of this, quoted in recent works, have more than 5 percent error. Treatises on cosmology 
deal with the uncertainty in the measured value of Hq by putting it as Hq = 100/z km s Mpc -1 and leaving h in the 
expressions as some undetermined parameter. The age of the universe has been another common stumbling block for 
the big bang theory. Many-a-time, it happens that the age of the universe in the big bang model is smaller than the ages 



1. INTRODUCTION 



of several celestial objects when measured directly. An easy way to solve this age problem is to include a small but 
positive cosmological constant A into the Einstein equation. The cosmological constant has a repulsive effect at large 
scales but at the same time it is assumed to have no other interactions with matter. (Energies with the latter property 
are sometimes termed as 'dark energy', when they are assumed to be present in substantial amount at present.) By its 
inclusion one obtains a larger age for the universe and thus can explain away the problem related to it. 

In spite of all these difficulties, many cosmologists believed that the big bang model provided a reliable and tested 
accounting of the universe from about 0.01 sec after the Bang until today, some 15 billion years later 1 1]. Thus in the 
1980s, cosmologists were primarily occupied with some conceptual issues with the big bang model, which are referred 
to as the 'cosmological problems'. The solution to some of these problems was found in the theory of inflation, an 
exponential expansion of the universe that might have occurred in its very early history, which can then lead to a flat 
universe. The energy which is supposed to have driven the inflationary expansion was a scalar field potential energy, 
which is akin to that of a cosmological constant, since it has an equation of state approximately of the form p\ = — Pa. 
But in the early 90s, there were no reasons to include a dark energy in the present universe. The euphoria regarding the 
validity of the big bang picture reached the pinnacle when the cosmic background explorer (COBE) satellite measured 
the fluctuations in the CMBR in 1992 to get its magnitude close to that predicted by the model. 

A recent development in theoretical cosmology is that the energy corresponding to the A-term or dark energy has 
again come back to the present universe scenario, as it appeared in the literature several times in the past for various 
reasons. Its present phase began with the measurement of the Hubble parameter reported in [2, 3], that gave a specific 
and high range h > 0.7. This, along with the flatness imposed by the theory of inflation, created a short spell of an age 
problem, which required a A-term for its solution. Then came the Type la Supernovae (SNe la) data |4, 5], which in 
turn required the presence of a A-term large enough to cause an accelerated expansion. The most recent release of SNe 
data | a Ells! prompts many cosmologists even to speculate that some extremely unphysical energy densities, such as 
those with w < — 1 in the equation of state p = wp are required to explain the data. 

Thus unlike in the past, it seems that the suspected dark energy in the universe cannot be regarded as a cosmological 
constant of general relativistic origin or as the vacuum energy encountered in quantum field theories. Our knowledge 
of the physical world then remains deficient since no tested theory involves such a dark energy. But in this talk I wish 
to caution that attempting to analyze the new SNe data by assuming exotic energy densities with strange equations 
of state itself is very misleading and argue that the reasonable remaining option is to analyze them cosmographically, 
without making any specific assumptions on the energy densities in the universe. Such an analysis of SNe data in |9] 
assumes only homogeneity and isotropy for the universe; i.e., the universe is assumed to have a Robertson- Walker 
(RW) metric. The scale factor of the universe can most naturally be expanded into a Taylor series in time about the 
present epoch and we attempt to find its coefficients from observation, as it was done to evaluate the Hubble parameter 
and deceleration parameter in the original cosmographic approach. Terms up to fifth-order were kept in the above series 
and it was assumed that they make a good approximation. The likelihoods for the various expansion rates obtained in 
our calculation are very broadly peaked and these indicate that there are a variety of choices for the energy densities 
in the present epoch. However, also this model-independent analysis confirmed that the universe is accelerating at 
present. 

An important question I try to answer in this talk is whether the data really endorse that the universe was decelerating 
in the past The answer has serious implications for modern theoretical cosmology for it is almost entirely built 
on the firm belief that the universe was decelerating in the past. Theories of the early universe, such as those of 
nucleosynthesis, CMBR, formation of structures, etc., were developed before the 1990s, under the assumption that the 
universe was always decelerating. Now with the arrival of the new SNe data, there appears to be a consensus that the 
universe is accelerating at present, and that this may be due to substantial amount of dark energy. But it should be 
noted that an accelerating universe was not in anyone's wildest imagination, at any time before these SNe observations 
were made. Obviously, this is sufficient reason to suspect the validity of the entire early universe theories. When the 
theoretical and observational status of the present universe has suffered such a setback, it would be naive to believe 
that the theories and observations regarding the early universe remain intact. The attempts of some cosmologists to 
deduce from the SNe data the redshift at which universe changed from deceleration to acceleration can only be viewed 
prima facie as a damage control measure. My attempt in this paper is to inspect this claim, whether the the latest SNe 
data give any indication that the universe changed from deceleration to acceleration at any foreseeable past! 

To answer our question, the Bayesian method of model comparison is adopted in which a model with some 
decelerating phase in the past is compared with another one having no such phase. In both cases, we assume the 
present universe to be accelerating. It is claimed that this method of comparison is more robust than other investigations 
seeking evidence for a decelerating past since both the Bayesian and cosmographic approaches mentioned above are 
used here. As suggested earlier, the scale factor is expanded into a fifth-order polynomial in time (fifth-order is required 



for sufficient accuracy) and then the combinations of various coefficients in this expansion were separated into those 
which correspond only to acceleration during the entire period and those which had at least some decelerating phase 
in this period. Considering these as rival models, the Bayes factor is calculated. The results show that if at all there 
is evidence, it is in favor of an always accelerating model; at best one can say that the SNe data are not capable of 
discriminating the "always accelerating" and "decelerating in the past" cosmological models. 

The talk is organized as follows. In Section 2, we state the standard methods to obtain the apparent magnitude- 
redshift relation of objects and the best fit parameter values in a model. Section 3 gives an overview of the Bayesian 
model comparison technique. It is in Section 4 that we present the Bayesian, model-independent analysis of our 
problem. The last section comprises a discussion of our results. 



2. THE APPARENT MAGNITUDE-REDSHIFT RELATION 

The conventional method of analyzing the apparent magnitude-redshift data relies on first solving the Einstein 
equations in cosmology for a{t), the scale factor appearing in the line element 
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Here p is the pressure and p is the density. We also need an equation of state of the form p = wp to solve the above. 
(w is called the equation of state parameter.) For each model, the solution of these equations helps us to predict the m-z 
relation as follows. Put ds = in Eq. Q for light coming from an SN, which is at position r\ at time t\, and reaches 
us at r = at the present time to and then integrate 
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This helps to find r\ or r\a.Q using the solution a(t) of Eqs. and obtained for the model under consideration. 
Then we compute the luminosity distance D = r\ «o ( 1 + z) and from it the apparent magnitude 
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Here we should make use of the relation 
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to obtain t\ in terms of z- It may be noted that D and hence m can be functions of M, «o, Hq, k, and other parameters in 
a(t) such as £l m = p m j p c , £2a = Pa/ Pc, etc., depending upon the model chosen. 

The evaluation of the best fit values of the parameters £2,„, £2a, etc., in a model is usually based on the minimum of 
the ^-statistic 
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where denotes the measured value of the apparent magnitude of the k SN, m/ ( (M,Ho,Q. m ,Q. 
values (from theory) and Cfy the uncertainties in the measurement of the observables. 
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3. BAYESIAN MODEL COMPARISON 



The Bayes's probability theorem tells us how to adjust our plausibility assumptions regarding a hypothesis when our 
state of knowledge changes through the acquisition of new data. Classical mathematicians such as Bernoulli, Bayes, 
Laplace and Gauss have found Bayes theorem useful in problems such as those in astronomy, thanks to its ability 
to learn. For example, Laplace has used this theory to estimate the masses of planets and to quantify the uncertainty 
of the masses due to observational errors. Later, because of the difficulties with assigning prior probabilities (which 
were mistakenly considered to be purely subjective expressions of a person's opinions about hypotheses), Bayes's 
probability theory has gone out of favor in physical sciences and was replaced by the more apparently objective 
'frequentist' probability approach. But the realization that the frequentist definition of probability is as subjective as 
the Bayesian has called forth a re-examination of this controversy. When compared to the improper application of 
frequentist probability theory, the Bayesian approach is most powerful in problems such as those in astronomy, where 
the process of learning and also the quantification and readjustment of plausibility by the scientific community are 
very important. 

In the Bayesian model comparison approach [11], we consider two models (say, M, and Mj) as rival hypotheses 
and try to compare them by evaluating odds ratios between the posterior (i.e., after analyzing the data) probabilities 
p(M | D,I) for these models, given the data D and also assuming the truth of some background information /. For this, 
the Bayes's theorem is made use of. This theorem states that the posterior probability for a hypothesis Hi given D and 
/is 

p^ D ^ = — — ■ (6) 

Here p(Hi\I) is called the prior probability for Hi, given the truth of / alone. p(D\Hi,I) is the probability for obtaining 
the data D if the hypothesis H, and / were true and is called the likelihood for the hypothesis. The factor in the 
denominator serves the purpose of normalization. 

Considering the truths of models as our hypotheses, the Bayes' theorem helps us to write the odds ratio between M, 
and Mj, which is the ratio between the corresponding posteriors, as 

p(Mj\D,I) = p(M,\I) P (D\M h I) ^ 
p(Mj\D,I) p(Mj\I)p(D\Mj,I) ,J - 

If the information / does not prefer one model over the other, the prior probabilities get canceled out and the odds 
ratio reduces to Bij, which is called the Bayes factor. 
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For parameterized models, the probability for the data given the truth of the model, also called the likelihood for the 
model, can be written as 



p{D\M h I)=^{Mi) = J da J dp...p{a,P,...\Mi)&i{a,p,...). 

Here, p(a,f5,...\Mi) is the prior probability for a . /3 . . . and the likelihood for the parameters in the model M,- is given 
by 

J^(a,j3,...) = exp[- Z , 2 (a,j3,..)/2]. (9) 

The commonly accepted interpretation of the Bayes factor is as follows: If 1 < B^ < 3, there is evidence against 
model Mj, but it is not worth more than a bare mention. If 3 < Bjj < 20, this evidence is positive. If 20 < Bij < 150, it 
is strong and if > 150, the evidence is very strong. 

As familiar examples of cosmological models, we consider the following for illustrating the Bayesian model 
comparison technique. 

( 1 ) Friedman-Lamaitre-Robertson- Walker (FLRW) model 

In this model, the universe contains two components for the total energy density p ; the first corresponds to ordinary 
matter and the second to a cosmological constant A. i.e., 



P = Pm + Pa, P = Pm + PA, 



A = constant, p\ = 



p m = 0, p A = -p A . 
(2) FLRW model with inflation (Flat universe) 

The same as the above with the additional condition £2 = £2,„ + £2 a = 1, which reduces the effective number of 
parameters by unity. 

Bayes factors between FLRW models with and without inflation (Models M\ and Mi respectively) obtained in full 
are 

B12 = 1 .6 using the data in 0] (Fit C) 

B12 = 1.1 using the 'All SCP SNe' data in 

This example shows that the inflationary model is slightly more favored than the general FLRW model, but the 
evidence is not worth more than a bare mention. When the Bayes factors are close to unity like this, it is highly 
objectionable to state that one or the other model is 'ruled out' by the data. It may also be reminded that the price for 
overlooking this fact would be frequent and unavoidable U-turns. 



4. THE MODEL-INDEPENDENT APPROACH TO THE PROBLEM 

In this section, we examine the claim that the universe changed from deceleration to acceleration at some time in the 
past, by adopting the Bayesian, model-independent approach. 

As stated in the introduction, the basic assumption in the model-independent or cosmographic approach is that the 
the scale factor of the universe can be approximated by the first six terms in a Taylor expansion I9 ul0ll . i.e., 
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We shall use this expression in the integral in Eq. (0). Parameters in this case are ao, h, qo, ro, so, uq, M, and k. 
We assume flat priors for the parameters in the ranges —19.6 < M < —19.1 magnitudes, 0.6 < h < 0.8, —2 < go < 0, 
-15 < r < 15, -65 < so < 65, and -150 < u < 150, 8000 > a > 3000 Mpc, and k = ±1. These ranges are as 
obtained in [9]. 

The two hypotheses we want to compare are (1) the universe is always accelerating from time to + T p to the present 
epoch fo (model M\) and (2) there is at least one decelerating epoch for the universe during this period (model M%). 
For any particular combination of parameter values, a sure test for the occurrence of deceleration during T p < T < 
is to plot 
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for this interval and to see whether it becomes negative at any time during the period. 

The values of the Bayes factor obtained while using the 54 "All SCP" SNe (data D\) in 0] (as reproduced in 01) 
and the 157 "gold" data points (data D2) in 1 8] are given in Table 1 . The results show that except in one case, the Bayes 
factors favor the 'always accelerating' Model 1. But in all cases, the evidences are weak. 



5. DISCUSSION 

We see that there is hardly any evidence in support of a decelerating past for the universe in the past 15 Gyrs, except in 
the case with T p = —2 x 10 17 s and data D^- These results are obtained using a combined Bayesian, model-independent 
approach. However, as per the interpretation of the Bayes factor given in Sec. 3, the evidences against a decelerating 
past for the universe are in general weak. Thus it may safely be concluded that these data alone cannot discriminate 
the two hypotheses we compared. In summary, we state that even the most precise cosmological data, which is that 



TABLE 1. Bayes factors between the "al- 
ways accelerating" (Mi) and "decelerating in 
the past" (M2) cosmological models. 
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of the SNe la, cannot provide any clue as to whether the universe changed from an earlier decelerating phase to the 
present accelerating one at some time in the past. 

Thus we are still at the starting point in cosmology, where the expansion of the universe is the only authentic piece 
of information, validated by known physics. 
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